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Abstract 

Measurements of inclusive W and Z boson production cross-sections in pp col- 
lisions at ^/s = 7TeV using final states containing muons are presented. The 
data sample corresponds to an integrated luminosity of 37pb~^ collected with 
the LHCb detector. The W and Z bosons are reconstructed from muons with a 
transverse momentum above 20GeV/c and pseudorapidity between 2.0 and 4.5, 
and, in the case of the Z cross-section, a dimuon invariant mass between 60 
and 120GeV/c2. The cross-sections are measured to be 831 ± 9 ± 27 it 29 pb for 
W+, 656 ± 8 ± 19 ± 23 pb for W and 76.7 ± 1.7 ± 3.3 ± 2.7 pb for Z, where 
the first uncertainty is statistical, the second is systematic and the third is due 
to the luminosity. Differential cross-sections, W and Z cross-section ratios and 
the lepton charge asymmetry are also measured in the same kinematic region. 
The ratios are determined to be (Tiy+_i.^+^,/cr^-_>^-p =1.27 it 0.02 it 0.01 and 
(c^w+^fi+u^ '^w-^tj.-iy)/'^z^ii^i =19.4 lb 0.5 it 0.9. The results are in general agree- 
ment with theoretical predictions, performed at next-to-next-to- leading order in 
QCD using recently calculated parton distribution functions. 
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1 Introduction 



The measurement of the production cross-sections for W and Z bosons constitutes an im- 
portant test of the Standard Model and provides valuable input to constrain the proton 
parton density functions (PDFs). Theoretical predictions are known to next-to-next-to- 
leading-order (NNLO) in perturbative quantum chromodynamics (pQCD). These calcula- 
tions are in good agreement with recent measurements at the LHC from the ATLAS [l|[2] , 
and the CMS [syi] experiments as well as with the results from the pp collider experi- 
ments at the SppS [5|[6] and the Tevatron [7 10 . The dominant theoretical uncertainty 



on the cross-sections arises from the present knowledge of the PDFs and the strong 
couphng constant. The accuracy strongly depends on the pseudorapiditjj^ range; con- 
sequently, measurements by LHCb, which is fully instrumented in the forward region 
2.0 < 1] < 5.0, can provide input to constrain the PDFs, both for pseudorapidities r] > 2.5 
and in the region which is common to ATLAS and CMS, 2.0 < rj < 2.5. Besides the 
determination of the W and Z boson cross-sections, the measurement of their ratios 

Rwz = {o'w+^tJ.+L' + o■w-^fJ.-p)/o^z^^J,^_l and Rw = CTw+^ti+v/'^w-^ii-u and of the W 
production charge asymmetry constitute important tests of the Standard Model, as ex- 
perimental and theoretical uncertainties partially cancel. The W charge asymmetry is 



sensitive to the valence quark distribution in the proton 11 and provides complementary 
information to the results from deep-inelastic scattering cross-sections at HERA [12] as 
those data do not strongly constrain the ratio of u over d quarks at low Bjorken x, where 
X is the proton momentum fraction carried by the quark. Measurements of W and Z 
boson production at LHCb have a sensitivity to values of x as low as 1.7 x 10^^ and 
will contribute significantly to the understanding of PDFs at low x and reasonably large 
four-momentum transfer Q^, which corresponds to the squared mass of the W or the Z 
boson. 

The measurements of the inclusive W and Z cross-section^ in pp collisions at a centre- 
of-mass energy of 7 TeV, using final states containing muons, are presented in this paper. 
The analysis is based on data taken by the LHCb experiment in 2010 with an integrated 
luminosity of 37pb~^. The cross-sections are measured in a fiducial region corresponding 
to the kinematic coverage of the LHCb detector, where the final state muons have a 
transverse momentum, p^^, exceeding 20GeV/c and lie within the pseudorapidity range 
2.0 < 7]^ < 4.5. This range is smaller than the LHCb acceptance in order to avoid 
edge effects for the acceptance. In addition, the invariant mass of the muons from the 
Z boson must be in the range 60 < M^^ < 120GeV/c^. Results are presented for the 
total cross-sections and cross-section ratios. Cross-sections are also measured in bins of 
muon pseudorapidity for W , and in bins of Z rapidity (y^) for Z production. Because 
of the presence of the neutrino, the production asymmetry between and W~ cannot 
be reconstructed as a function of the boson rapidity. Instead it is measured as a function 
of the experimentally accessible muon pseudorapidity, rj^, and referred to as the lepton 



"'^The pseudorapidity rj is defined to be 77 = — lntan(0/2), where the polar angle is measured with 
respect to the beam axis. 

^Throughout this paper Z includes both the Z and the virtual photon (7"^) contribution. 
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charge asymmetry = {aw+^fi+u - (Jw-^iM-u)/{(^w+^fM+u + crw-^^i-u)- To constrain 
the PDFs, it is useful to measure A for different pi^ thresholds. The data are compared 
to NNLO and NLO pQCD predictions with recent parametrisations for the PDFs. The 
signal efficiency and background contribution are mostly derived from data. 

The remainder of the paper is organised as follows. Section 2 describes the LHCb 
detector and the Monte Carlo samples. Section 3 describes the selection of the W and Z 
candidates, the backgrounds, the determination of the purity and the signal efficiencies. 
The measurement of the cross-sections as well as the systematic uncertainties are discussed 
in Sect. 4. The results are presented in Sect. 5 and conclusions in Sect. 6. 



2 LHCb detector and Monte Carlo samples 

The LHCb detector is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < 77 < 5, designed for the study of particles containing 6 or c quarks. The detec- 
tor includes a high precision tracking system consisting of a silicon-strip vertex detector 
(VELO) surrounding the pp interaction region, a large-area silicon-strip detector (TT) 
located upstream of a dipole magnet with a bending power of about 4Tm, and three 
stations of silicon-strip detectors (IT) and straw drift-tubes (OT) placed downstream. 
The combined tracking system has a momentum resolution Ap/p that varies from 0.4% 
at 5GeV/c to 0.6% at lOOGeV/c, and an impact parameter (IP) resolution of 20 |xm for 
tracks with high transverse momentum. Charged hadrons are identified using two ring- 
imaging Cherenkov detectors. Photon, electron and hadron candidates are identified by 
a calorimeter system consisting of scintillating-pad and pre-shower detectors, an electro- 
magnetic calorimeter and a hadronic calorimeter. Muons are identified by a muon system 
composed of alternating layers of iron and multiwire proportional chambers. The trigger 
consists of a hardware stage, based on information from the calorimeter and muon sys- 
tems, followed by a software stage which applies a full event reconstruction. To avoid the 
possibility that a few events with high occupancy dominate the CPU time of the software 
trigger, a set of global event cuts is applied on the hit multiplicities of most subdetectors 
used in the pattern recognition algorithms. 

Several Monte Carlo (MC) simulated samples are used to develop the event selection, 
estimate the backgrounds, cross-check the efficiencies and to account for the effect of the 



underlying event. The Pythia 6.4 14 generator, configured as described in Ref. 15 
with the CTEQ611 16 parametrisation for the PDFs is used to simulate the processes 
Z — 7- /i/i, Z —7- rr, W ^ fiu and W — > tu. The hard partonic interaction is calculated 
in leading order pQCD and higher order QCD radiation is modelled using initial and 



final state parton showers in the leading log approximation 17 . The fragmentation into 



hadrons is simulated in Pythia by the Lund string model 18 . All generated events are 



passed through a Geant4 19 based detector simulation, the trigger emulation and the 
event reconstruction chain of the LHCb experiment. 

Samples of W ^ and Z /i/i simulated events with one muon in the LHCb 
acceptance have been reweighted to reproduce the NNLO pi^ distribution. These samples 
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are referred to as IV-MC and Z-MC, respectively. In the first step a correction factor 
is calculated as a function of the generated muon transverse momentum by determining 
the ratio of the generated spectrum, as simulated by the Powheg 20 generator at 



NLO, to the generated pf^ spectrum from Pythia. In the second step the events are 
reweighted with a factor given by the ratio between the NNLO and NLO prediction as 
calculated with Dynnlo 21 . This factor is calculated as a function of the rapidity of 



the boson. As an alternative, Pythia samples have been reweighted to reproduce the pi^ 
distribution as calculated with Resbos 22 . Resbos includes a NLO calculation plus 



next-to-next-to-leading-log resummation of QCD effects at low transverse momentum. 



3 Selection of W and Z events 

3.1 Muon reconstruction and identification 

Events with high transverse momentum muons are selected using a single muon trigger 
with a threshold of pf^ > lOGeV/c. Tracks are reconstructed starting from the VELO, 
within which particle trajectories are approximately straight, since the detector is located 
upstream of the magnet. Candidate tracks are extrapolated to the other side of the 
magnet and a search is made for compatible hits in the IT and OT sub-detectors. An 
alternative strategy searches for track segments in both the VELO and IT/OT detectors 
and extrapolates each to the bending plane of the magnet, where they are matched. 
Once VELO and IT/OT hits have been combined, an estimate of the track momentum 
is available and the full trajectory can be defined. Finally, hits in the TT sub-detector 
are added if consistent with the candidate tracks. Thus, the presence of TT hits can be 
considered as an independent confirmation of the validity of the track. 

Muons are identified by extrapolating the tracks and searching for compatible hits in 
the four most downstream muon stations. For the high momentum muons that concern 
this analysis, hits must be found in all four muon stations. In total, the muon candidate 
must have passed through over 20 hadronic interaction lengths of material. 



3.2 Selection of Z — )- candidates 

Candidate Z ^ fifi events are selected by requiring a pair of well reconstructed 
tracks identified as muons; the invariant mass of the two muons must be in the range 
60 < M^^ < 120 GeV/c^. Each muon track must have p!^ > 20 GeV/c and lie in the range 
2.0 < 7]^ < 4.5. The relative uncertainty on the momentum measurement is required to 
be less than 10% and the probability for the x^/ndf for the track fit larger than 0.1%, 
where ndf is the number of degrees of freedom. In total, 1966 Z candidates are selected; 
their mass distribution is shown in Fig. [TJ The data are not corrected for initial or final 



state radiation. A Crystal Ball 23 function for the Z peak, and an exponential distri- 
bution for both the off-resonance Drell-Yan (7*) production and the small background 
contribution are fitted to the distribution. The fitted mass 90.7 ± 0.1 GeV/c^ and width 
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3.0 ± O.lGeV/c^, where the uncertainties are statistical, are consistent with expectation 
from simulation. 




Figure 1: Invariant mass of the selected muon pairs. The fitted distribution to the data is 
shown as a solid line and the contribution from background and off-resonance Drell-Yan 
production dashed line. 



3.3 Z — )- event yield 

The background contribution to the Z — /x/x analysis is very low. Five different 
sources are investigated. 

1. Decays from Z ^ tt contribute, if both taus decay leptonically to muons and 
neutrinos. The tau background is estimated from simulation, with the Z cross- 
section fixed to the cross-section measured in this analysis, to contribute 0.6 ±0.1 
events to the total sample. 

2. Decays of heavy flavour hadrons contribute to the background if they decay semilep- 
tonically ("heavy flavour" background). The contribution is estimated from a sam- 
ple, which is enriched in background. "Non-isolated" muons are selected with 
•p^ > 15 GeV/ c and M^^ > 40 GeV/ and the scalar sum of the transverse mo- 
menta of all tracks in a cone of half angle 0.5 in ■)]-(/) around the muons larger than 
4GeV/c; here is the azimuthal angle measured in radians. A fit to the invariant 
mass distribution at low masses is then used to estimate the background contri- 
bution in the Z mass region. The heavy flavour contribution is estimated to be 
3.5 ± 0.8 events. 
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3. Pions or kaons may be misidentified as muons if they decay in flight ("decay-in- 
flight" background) or if they travel through the calorimeters and are identified 
by the muon chambers ("punch-through" background). This background should 
contribute equally in same-sign and opposite-sign combinations of the muon pair. 
No event is found in the Z selection with both tracks having the same charge. The 
contribution from muon misidentification is estimated to be less than one event. 

4. W pair production contributes to the sample if both W bosons decay to a muon 
and a neutrino. This contribution corresponds to 0.2 ±0.1 events as estimated with 
Pythia MC simulation. 

5. Decays of top quark pairs may contribute if both top quarks decay semileptonically. 
Pythia MC simulation predicts a contribution of 0.5 ± 0.2 events. 

The total background contribution in the Z sample in the range 60-120 GcV/c^ 
amounts to 4.8 ± 1.0 events. This corresponds to a purity = 0.997 ± 0.001. The 
purity is defined as the ratio of signal to candidate events. No significant dependence on 
the boson rapidity is observed. 

3.4 Selection of W ^ iiu candidates 

In leading order QCD, W ^ fiv events are characterised by a single high transverse 
momentum muon that is not associated with other activity in the event. As only the muon 
can be reconstructed in LHCb, the background contribution is larger for the W than for 
the Z candidates. Therefore, more stringent requirements are placed on the track quality 
of the muon and additional criteria are imposed in order to select W candidates. 

The optimisation of the W selection and the evaluation of the selection efficiency 
make use of a "pseudo-VF" control sample obtained from the previously described Z 
selection, where each of the muons is masked in turn, in order to mimic the presence of 
a neutrino and fake a, W ^ fiu decay. Excellent agreement is observed for all variables 
of interest between pseudo-W^ and W simulated samples with the exception of those that 
have an explicit dependence on the transverse momentum of the muon, as the underlying 
momentum distribution differs for muons from Z and W. 

The identification oiW ^ nu candidate events starts by requiring a well reconstructed 
track which is identified as a muon. The track must have a transverse momentum in the 
range 20 < < 70GeV/c within a pseudorapidity range 2.0 < 1]^^ < 4.5. The relative 
error on the momentum measurement must be less than 10%, the probability for the 
X^/ndf of the track fit must be greater than 1%, and there must be TT hits associated 
to the track. The last requirement reduces the number of combinations of VELO and 
IT/OT information that have been incorrectly combined to form tracks. 

To suppress background from Z — )■ pp decays, it is required that any other identified 
muon in the event has a transverse momentum below 2GeV/c. This removes the events 
where both muons have entered the LHCb acceptance. 
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Figure 2: Distributions for p™"^" (top) and E^^'^ (bottom). The points are for muons 
from pseudo-H^ data, the yellow (shaded) histograms are for VT-MC simulation, while the 
open histograms are for muons from QCD background with IP > 80 (xm from data. All 
distributions are normalised to unity. 



Identified muons can originate from background processes of heavy flavour decays, 
or misidentification of pions and kaons due to decay-in-fiight or punch-through ("QCD 
background"). In all such cases, the identified muon is usually produced in the same 
direction as the other fragmentation products, in contrast to muons from W decays which 
tend to be isolated. The isolation of the muon is described using the charged transverse 
momentum, and neutral transverse energy, E^"^^, in a cone around the candidate 
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muon. The quantity p^'^^ is defined as the scalar sum of the transverse momentum 
of all tracks, excluding the candidate muon, satisfying iy(A0p~+~(Ar7^ < 0.5, where 
A0 and Arj'^ are the differences in and t] between the muon candidate and the track. 
The quantity E^^'^ is defined in a similar way, but summing the transverse energy of 
all electromagnetic calorimeter deposits not associated with tracks. The distributions 
for j9™"^ and E^'^^ are shown in Fig. [2] for pseudo- W data, VT-MC and muons with 
pi^ > 20 GeV/ c and an IP larger than 80 \xm. The IP of the muon is defined as the distance 
of closest approach to the primary vertex calculated from the other tracks in the event 
excluding the muon candidate. The sample with high IP is enriched with muons from 
decays of heavy flavour hadrons, showing the typical shape of QCD background. There 
is agreement between pseudo- W data and H^-MC, while the shape for the heavy flavour 
events is quite different. To suppress QCD background, it is required that pJj?"*^ < 2 GeV/c 
and < 2 GeV. 

Muons originating from semi-leptonic decays of heavy flavour hadrons can be further 
suppressed by a cut on the IP. Due to the lifetimes of the B and D mesons, these muons 
do not originate from the primary pp interaction. The IP distribution is shown in Fig. [3] 
for pseudo- W events, W-MC, and simulated semi-leptonic decays of hadrons containing a 
6 or c quark. The pseudo- W events and VT-MC are in agreement and peak at low values 
of IP, in contrast to the heavy flavour background. For the W candidate selection it is 
required that IP < 40 [im. This cut also removes a large fraction of the background from 
W ^ TV and Z ^ tt decays. 

Pions and kaons that punch-through to the muon chambers can be distinguished 
from true muons as they leave substantial energy deposits in the calorimeters. Figure |4] 
shows the summed energy, E, in the electromagnetic and hadronic calorimeter associated 
with the particle, divided by the track momentum, p, for pseudo-Ty events, M^-MC, and 
hadrons with > 20GeV/c in randomly triggered events. By requiring E/pc < 0.04 the 
punch-through contamination can be reduced to a negligible level. The disagreement be- 
tween pseudo-py data and simulated VT-MC in Fig. |4]is caused by the different underlying 
momentum distribution for muons from W and Z. 

3.5 W fjLiy event yield 

After the W selection requirements are imposed 14 660 and 11 618 W~ candidate 
events are observed. The W ^ fiu signal yield has been determined by fitting the pf^ 
spectra of positive and negative muons in data, to template shapes for signal and back- 
grounds in five bins of rj'^. The fit is performed with the following sources for signal and 
background with the shapes and normalisations as described below. 

1. The W ^ fiu signal template is obtained using the VT-MC. The normalisation is 
left free to vary in each bin of rj'^ and for each charge. 

2. The shape of the template of the largest background, Z — )■ fifi, is taken from the 
Z-MC. The normalisation is fixed from data by counting the number of Z events, 
scaled by the ratio of events with one muon in the LHCb acceptance to events with 
both muons in the acceptance, as determined from Z-MC. The ratio is corrected for 
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Figure 3: Muon IP distribution for pseudo-VK events as points, ly-MC as a yellow 
(shaded) histogram, and muons from simulated semi-leptonic decays of hadrons contain- 
ing a b quark in the full open histogram or a c quark in the dashed open histogram. All 
distributions are normalised to unity. 




Figure 4: E/pc for pseudo-VK events as points, VT-MC as a yellow (shaded) histogram, 
and for hadrons from randomly triggered events in the open histogram. The energy E is 
the sum of the energies in the electromagnetic and hadronic calorimeter associated with 
the particle. All distributions are normalised to unity. 
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the different reconstruction and selection efficiencies for W and Z as derived from 
data. This gives an expectation of 2435 ± 101 background events ((9.3 ± 0.4)% of 
the total sample) in good agreement with 2335 ± 25 events found from simulation. 

3. The shape of the W —> tu and Z tt templates are taken from Pythia. The 
Z — )■ rr template is scaled according to the observed number of Z events. These r 
backgrounds constitute 2.7% of the total sample. 

4. The heavy flavour template is obtained from data by requiring that the muon is 
not consistent with originating from the primary vertex (IP > 80|J.m). The nor- 
malisation is determined from data applying all requirements except for the impact 
parameter and fitting the resulting IP distribution to the two templates shown in 
Fig. [3] the pseudo-W^ data to describe the signal, and the simulated heavy fiavour 
events to describe the background. The heavy fiavour contribution is estimated to 
be (0.4 ± 0.2)% of the total sample. 

5. The punch-through contribution from kaons and pions is largely suppressed by the 
requirement on E/pc. The E/pc distribution in Fig.|4]is fitted to pseudo-VT data for 
the signal, and a Gaussian for the punch-through, in order to estimate the punch- 
through contribution. This is found to be negligible (0.02 ± 0.01)% of the total 
sample, and also has a shape very similar to the decay-in-fiight component. Hence, 
this component is not considered when determining the signal yield. 

6. The decay-in-fiight shape is found from data in a two-step procedure using all events 
selected throughout 2010 by any trigger requirement. First, tracks with a transverse 
momentum between 20 and 70 GeV/ c are taken to describe the p^ spectrum of 
hadrons; tracks that fired a muon trigger are excluded from the sample. Second, 
this spectrum is weighted according to the probability for a hadron to decay-in- 
fiight. This probability is defined as the fraction of tracks identified as muons in 
randomly triggered events and is parametrised as a function of the momentum, p, 
by a function of the form 

1 - e-"/P, (1) 

as would be expected for a particle whose mean lifetime in the laboratory frame 
scales with its boost. Consistent values for a are found in each pseudorapidity 
bin and are in agreement with a calculation of the decay probability based on the 
mean lifetimes for charged pions and kaons, and the distance to the electromagnetic 
calorimeter before which the hadron must have decayed. The average of the deter- 
minations in each pseudorapidity bin defines the central value for a. The relative 
normalisation of positively to negatively charged tracks in each bin of pseudorapidity 
is fixed to that observed in randomly triggered events, but the overall normalisation 
in each bin of pseudorapidity is left free. 

The default fit has 15 free parameters: five parameters for the normalisation of in 
each of the pseudorapidity bins, five parameters for W~, and five parameters for the 
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Figure 5: Distribution of muon for positively (left panel) and negatively (right panel) 
charged muons in W candidate events, for the total fiducial cross-section (a). The plots 
(b) to (f) give the same information for the different rj'^ bins. The data (points) are 
compared to the fitted contributions from W~ and (light shaded). The background 
contributions are, from top to bottom in the legend: decay- in- flight, Z — )■ /x/i, r decays 
of W and Z, and heavy flavour decays. 
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contribution coming from the decay- in-flight. The normahsation of the other sources is 
fixed. The result of the fit is shown in Fig. [s] Integrated over both charges and pi^ it is 
found that (44.3 ± 1.2)% of the total sample is due to W+, (34.9 ± 1.1)% due to W, 
(8.5 ± 0.8)% due to the decay-in-flight contribution and the remainder due to the other 
backgrounds. The x^/^df of the fit is 1.002. The fit is repeated with the Z — )■ /i/i and 
W fxu template corrected with Resbos instead of NNLO, yielding (43.6 ± 1.2)% for 
W+ and (34.4 ± 1.1)% for W' with xV^df = 0.983. The average of the two fits, which 
gives a purity p^'' = 0.788 ± 0.021 for W+ and p^' = 0.784 ± 0.025 for W , is taken for 
the final result; half of the difference is taken as the systematic uncertainty. 



4 Cross-section measurement 



4.1 Cross-section definition 



Cross-sections are quoted in the kinematical range defined by the measurements. The 
cross-sections are measured in bins of r)'^ for the W and in bins of in case of the Z. 
The cross- section in a given bin of [t]^) is defined as 



crz- 



where N^{r]'^,r]j) is the number of Z candidates in the respective y^ bin with the two 
muons in the bins t]^ and rjj being reconstructed with the efficiency B^{ri^, rjj). Similarly, 
is the number of W candidates with the muon in the t]'^ bin. The purity of the sample 
(p^^^^), the acceptance (A^^^^), the correction factor for final state radiation (FSR) 
ifpsR^) the efficiency (e^) are determined per bin; C is the integrated luminosity. 
The total cross-section is obtained by summing the contributions of the five y^ or t]'^ bins. 



4.2 Signal efficiencies 

The data are corrected for efficiency losses due to track reconstruction, muon identifica- 
tion, and trigger requirements for both analyses. There is an additional selection efficiency 
in the W analysis due to the requirements on the number of additional muons, IP, E/pc, 
p^^^, E^'^'^, and on TT hits. All efficiencies are determined from data. 

The efficiencies for track reconstruction and muon identification are obtained using a 
tag-and-probe method in the Z sample. One of the muons in the Z sample (tag) satisfies 
all the track criteria. The other muon (probe) is selected with looser criteria that depend 
on the efficiency to be measured. The invariant mass of the dimuon candidates, recon- 
structed from the tag and the probe muons, must lie in the window of 20 GeV/c^ around 
the nominal Z mass. The tracking efficiency, which accounts for track reconstruction and 
the track quality requirements, is studied using well reconstructed tracks in the muon 
stations which are linked to hits in TT. The average track finding efficiency is about 90% 
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in the Z sample and about 86% for the muon in W events. The tracking efficiency for 
W is lower due to the more restrictive cuts on the track quality. The muon identification 
efficiency is determined with tracks without the muon identification requirement for the 
probe muon. The average single muon efficiency is above 99%. Both the tracking and the 
muon identification efficiencies agree with simulation within errors. 

The trigger efficiency contains two components, the ffist due to the efficiency of the 
single muon trigger and the other due to the global event cuts (GEC). The single muon 
trigger efficiency is determined using the Z — )■ /iyU sample. One muon is required to &e 
the single muon trigger. The trigger response of the other muon then defines the trigger 
efficiency. The requirement on the occupancy of the events depends on the multiplicity 
of the primary interactions. It was checked with a sample which did not have the GEC 
applied, that no events are lost if there is only one primary vertex reconstructed. The 
GEC efficiency as a function of the number of primary vertices is determined by adding 
randomly triggered events to events with Z iyV^ candidates with one primary vertex; 
on average it amounts to 93%. The overall trigger efficiency is calculated for each event 
depending on the lepton pseudorapidity and the primary vertex multiplicity. It is found 
to be about 88% for the Z and 75% for the W sample. 

The W selection requires that there are no other muons with pt^ > 2 GeV/ c, 
p^"""" < 2 GeV/c, E^p"'' < 2 GeV, IP < 40 [xm, and E/pc < 0.04. The selection efficiency is 
determined from the fraction of pseudo- W events that pass these requirements. A similar 
method is used to evaluate the efficiency for the requirement of TT hits associated to the 
muon track of the W candidate. Simulation studies show that with the exception of the 
E/pc distribution, the pseudo- data provide a consistent description of W fiu sim- 
ulation, as shown in Figs. [2] and |3j However, the harder muon spectrum in pseudo- W 
data leads to slightly lower values of E/pc than for muons produced in W decays. The 
simulation is used to determine this difference, which is only significant for rj'^ between 
2.0 and 2.5, where the efficiency for W events is estimated to be 2.1% lower than for 
pseudo- ly data. The selection efficiency is about 67% for 2.5 < t]'^ < 4.0 and drops to 
about 53% and 33% for the two bins at the edge of the acceptance with 2.0 < r/^ < 2.5 
and 4.0 < rj^ < 4.5, respectively. 

All the efficiencies have been checked for possible dependences on pf^, the azimuthal 
angle of the muon, magnet polarity, and r]'^. Only the latter exhibits a significant depen- 
dence, which is taken into account. Since any charge bias of the efficiencies would directly 
influence the measurement of the lepton charge asymmetry, it was checked there is no 
significant charge dependence within the uncertainties of the efficiencies. The efficiency 
corrections are applied as a function of the pseudorapidity of the muons except the GEC. 

The efficiencies are un correlated between pseudorapidity bins but correlated for W^, 
W~ and Z. These correlations are taken into account for the measurement of the lepton 
charge asymmetry and the cross-section ratios. 
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4.3 Acceptance 



The selection criteria for the W and Z define the fiducial region of the measurement. 
Simulated events are used to determine the acceptance A, defined as ^ = Nj-ec/Ngen- 
Here, Nj-ec is the number of reconstructed events satisfying the cuts on the pseudorapidity 
and the minimal momenta of the reconstructed muons, as well as on the dimuon mass in 
the case of the Z analysis. Similarly, A^^gen is the number of generated events with the cuts 
applied on the generated muons. The acceptance is estimated with ly-MC and Z-MC. It 
is found to be consistent with unity for the Z and above 0.99 for the W analysis. For the 
latter, the acceptance corrects for the small loss of events with pi^ > 70 GeV/ c. 



4.4 Luminosity 

The absolute luminosity scale was measured at specific periods during the data taking. 



using both Van der Meer scans 24 where colliding beams are moved transversely across 



each other to determine the beam profile, and a beam-gas imaging method 25 ,26 . For the 
latter, reconstructed beam-gas interaction vertices near the beam crossing point determine 
the beam profile. Both methods give similar results and are estimated to have a precision 
of order 3.5%. The knowledge of the absolute luminosity scale is used to calibrate the 
number of tracks in the VELO, which is found to be stable throughout the data-taking 
period and can therefore be used to monitor the instantaneous luminosity of the entire 
data sample. The dataset for this analysis corresponds to an integrated luminosity of 
37.1 ± 1.3pb-\ 



4.5 Corrections to the data 

The measured cross-sections are corrected to Born level in quantum electrodynamics 
(QED) in order to provide a consistent comparison with NLO and NNLO QCD pre- 
dictions, which do not include the effects of FSR. Corrections have been estimated using 



Photos 27 interfaced to Pythia. The Pythia spectrum of the electroweak boson 



has been reweighted to the NNLO spectrum as determined with Dynnlo 21 . The cor- 
rection is taken as the number of events within the fiducial cuts of the measurements after 
FSR divided by the number of events generated within the fiducial cuts. 

Pythia simulation is used to study bin-to-bin migrations for r]'^ and . No significant 
net migration is observed and no correction is applied. 



4.6 Systematic uncertainties 

Aside from the uncertainty on the luminosity measurement, the main sources of experi- 
mental uncertainties come from the efficiency determinations and the background estimate 
in the W analysis. The following sources have been considered: 

1. The relative uncertainties of the tracking, muon identification, trigger and GEC 
efficiencies are added in quadrature. They lead to a systematic uncertainty for the 
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Table 1: Contributions to the systematic uncertainty for the total Z and W cross-sections. 
The different contributions are discussed in Sect. 14.61 



Source 


Aaz^^^ (%) 


AcTHZ+^/t+i/ (%) 


Aavy-^^-p (%) 


Signal purity 


±0.1 


±1.2 


±0.9 


Template shape (fit) 




±0.9 


±1.0 


Efficiency (trigger, tracking, muon id) 


±4.3 


±2.2 


±2.0 


Additional selection 




±1.8 


±1.7 


FSR correction 


±0.02 


±0.01 


±0.02 


Total 


±4.3 


±3.2 


±2.9 


Luminosity 


±3.5 


±3.5 


±3.5 



total cross-sections of 4.3% (2.2%, 2.0%) for the Z (iy+, W). 

The statistical uncertainty on the efficiency of the additional selection cuts for the 
W analysis translates into a 1.8% (1.7%) systematic uncertainty on the total 
(iy~) cross-section. 

The uncertainty of the background contribution for the Z analysis is small; the 
uncertainty in the determination of the sample purity leads to a 0.1% uncertainty 
on the total cross-section. 

Both the shape and normalisations of the templates used in the W fit are considered 
as an additional source of uncertainty. To determine this systematic uncertainty 
each of the following sources is varied in turn, the data are refitted to determine the 
fraction of and W~ events, and the deviations from the original signal yield are 
combined in quadrature. The following variations are made: 



the difference of the two fits using different W and Z templates (see Sect. 3.1) 
leads to a variation on the W fractions of 0.8%; 

the normalisation of the Z component was changed by the statistical uncer- 
tainty with which it was determined, leading to a variation in the W fractions 
of 0.3%; 

the normalisation of the W ^ tv template was changed by the statistical 
uncertainty with which it was determined, leading to a negligible change in 
the W fractions, since this template shape is very similar to the decay-in-flight 
template which is allowed to vary in the fit; 

the heavy flavour template has also been changed by the statistical uncertainty 
with which it was determined leading to a negligible change in the W fractions; 

instead of leaving the relative normalisation of the decays-in-flight template 
between pseudorapidity bins to be free in the fit, this is fixed to the values 
observed in randomly triggered events, and the full fit performed with a single 
free parameter for the background; the W fractions change by 0.2%; 
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• the shape of the decay-in-flight template has been changed using different val- 
ues for a (see Eq. [T]) to describe the decay probability, corresponding to differ- 
ent regions in which the hadron must have decayed]^ no difference in the W 
fractions is observed. 

5. The uncertainty on the FSR correction is evaluated for each bin as the maximum 
of the statistical uncertainty of the correction factor and the difference between the 
weighted and unweighted FSR correction factor. 

The sources of systematic uncertainties are summarised in Table [T| together with the size 
of the resultant uncertainty on the W and Z total cross-sections. The total systematic 
uncertainty is the sum of all contributions added in quadrature. 



5 Results 

The inclusive cross-sections for Z — )■ /x/i and W ^ [iv production for muons with pl^ > 
20 GeV/c in the pseudorapidity region 2.0 < r)^ < 4.5 and, in the case of Z, the invariant 
mass range 60 < M^^ < 120 GeV/ are measured to be 

az^fif, = 76.7 ± 1.7 ± 3.3 ±2.7pb 
aw+^ti+u = 831±9±27±29pb 
c^VF-->M-p = 656±8±19±23pb, 

where the first uncertainty is statistical, the second systematic and the third is due to the 
luminosity. All the measurements are dominated by the luminosity and the systematic 
uncertainty. The latter is dominated by the limited number of events for the background 
templates and in the determination of the efficiencies. 

The ratios Rw = aw+^,i+u/o-w-^ti-p and Rwz = {aw+^ii+u + o-w-^,i-u)/crz^t^/i are 
measured to be 

Rw = 1.27 ±0.02 ±0.01 
Rwz = 19.4 ±0.5 ±0.9. 

Here, the uncertainty from the luminosity completely cancels. The systematic uncer- 
tainty on the trigger, muon identification, tracking and selection efficiencies, as well as 
the uncertainty on the purity are assumed to be fully correlated between and W~ . 
No correlation is assumed between the rj'^ bins, except for the purity. The uncertainty 
on the Z cross-section from the reconstruction efficiency is correlated between boson ra- 
pidity bins. The correlation of the uncertainty on the efficiencies between W and Z are 
estimated with MC simulation to be 90%. The full correlation matrix is given in the 
Appendix (Table [2]). The ratio of the W to Z cross-section is measured, for each charge 
separately, to be 

^Three different decay regions have been considered: from the interaction point, from the VELO and 
from the TT stations up to the electromagnetic calorimeter. 
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Figure 6: Measurements of the Z, W'^ and W~ cross-section and ratios, data are shown as 
bands which the statistical (dark shaded/orange) and total (light hatched/yellow) errors. 
The measurements are compared to NNLO and NLO predictions with different PDF sets 
for the proton, shown as points with error bars. The PDF uncertainty, evaluated at the 
68% confidence level, and the theoretical uncertainties are added in quadrature to obtain 
the uncertainties of the predictions. 



(Jw+^^l+v/(^z^^Ji^, = 10.8 ± 0.3 ± 0.5 
aw-^^^-9/(Tz^^,f, = 8.5 ±0.2 ±0.4. 

A summary of the measurements of the inclusive cross-sections aw+^fi+u, <^w-^^i-u 
and az^^i^i, and the ratios is shown in Fig. |6} The measurements are shown as a band 
which represents the total and statistical uncertainties. 

The results are compared to theoretical predictions calculated at NNLO with the 
program Dynnlo |21] for the NNLO PDF sets of MSTW08 [28], ABKM09 [29], JR09 |30), 
HERA15 (12) and NNPDF21 |3l) and at NLO for the NLO PDF set CTEQ6m [16] The 



scale uncertainties are estimated by varying the renormalisation and factorisation scales by 



' Dynnlo sets as to the value of ag at the mass of the Z boson as given by the different PDF sets. 
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Figure 9: Lepton charge asymmetry = {aw+^f^+u - <yw-^fi'u) / {crw+^^,+u + crw-^^i-u) 
in bins of muon pseudorapidity. The dark shaded (orange) bands correspond to the sta- 
tistical uncertainties, the hght hatched (yellow) band to the statistical and systematic 
uncertainties added in quadrature. Superimposed are NNLO (NLO) predictions as de- 
scribed in Fig [7j The MSTW08 values for r]^ < 2 represent the central value of the 
prediction. 



factors of two around the nominal value, which is set to the boson mass. The uncertainties 
for each set correspond to the PDF uncertainties at 68% and the scale uncertainties added 
in quadrature]^ 

While the W~ and Z cross-sections are well described by all predictions, the 
cross-section is slightly overestimated by the ABKM09 and NNPDF21 PDF sets. The ra- 
tio of the W~ to Z cross-sections agrees reasonably well with the predictions, but the 
to the Z ratio is overestimated by most of the predictions. The systematic uncertainties 
for the Rw almost cancel and also the theoretical uncertainties are much reduced. The 
Rw measurement tests the Standard Model predictions with a precision of 1.7% which 
is comparable to the uncertainty of the theoretical prediction. The ABKM09 predic- 
tion overestimates this ratio while all the other predictions agree with the measurement. 
Differential distributions are measured in five bins in for the Z and of r}^ for the W. 
Figure [7] shows the differential cross-section as a function of the rapidity of the Z boson 
together with NNLO (NLO) predictions with different parametrisation for the PDFs of 
the proton. The predictions agree with the measurements within uncertainties though 
all the predictions are lower than the measured cross-section for 2.5 < t]^ < 3.0. The 
differential cross-sections are listed in Table [3] in the Appendix. 

^The uncertainties for the PDF set from CTEQ6m which is given at 90% CL are divided by 1.645. 
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Figure 10: Rw = (yw+^^i+v/(^w-^^i'u in bins of muon pseudorapidity. The dark shaded 
(orange) bands correspond to the statistical uncertainties, the hght hatched (yellow) band 
to the statistical and systematic uncertainties added in quadrature. Superimposed are 
NNLO (NLO) predictions with different parametrisations as described in Fig [T] The 
MSTW08 values for r/^ < 2 represent the central value of the prediction. 



The differential distribution of the and W~ cross-section, the lepton charge asym- 
metry and the ratio Rw as a function of the muon pseudorapidity are shown in 
Figs, [sj |9] and 10 and listed in Tables |4] to [6] as a function of p'^. The measurement of 
the charge asymmetry and the W ratio provides important additional information on the 



PDFs particularly on the valence quark distributions 11 



Since the inclusive cross-section for is larger than for W~ , due to the excess of u 
over d quarks in the proton, the overall asymmetry is positive. The asymmetry and the 
W cross-sections strongly vary as a function of the pseudorapidity of the charged lepton, 
and even changes sign, owing to differing helicity dependence of the lepton couplings 
to the boson. This behaviour is reflected in the differential W cross- sections, where at 
large muon pseudorapidities the W~ cross-section is higher than the cross-section, as a 
consequence of the V — A structure of the W to lepton coupling. The cross-section and the 
asymmetry measurements are compared to the NNLO (NLO) predictions with different 
parameterisation for the PDFs. The ABKM09 prediction overestimates the measured 
asymmetry in three of the five bins. The other predictions describe the measurement 
within uncertainties. 

The asymmetry is also measured for two higher p'^ thresholds for the muons, at 25 
and 30 GeV/c. The result is shown in Fig. 11 and listed in Table |5} The NNLO prediction 
with MSTW08 parametrisation for the PDF also describes the measured asymmetry with 
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Fi gure 11: Lepton charge asymmetry A.^ — (^(J\y+_^^+i, — crw^^fi^u) / {'^w+-¥fj.+u~^(^w^-¥fj.~'u) 
for muons with pt^ >25 (top) and 30GeV/c (bottom), respectively in bins of muon pseu- 
dorapidity. The dark shaded (orange) bands correspond to the statistical uncertainties, 
the light hatched (yellow) band to the statistical and systematic uncertainties added in 
quadrature. The statistical uncertainty is undistinguishable from the total uncertainty. 
Superimposed are the NNLO predictions with the MSTW08 parametrisation for the PDF. 
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the higher cuts on the transverse momentum of the muon. 



6 Conclusions 

Measurements of inclusive W and Z boson production in pp coUisions at ^/s — 7 TeV with 
final states containing muons have been performed using 37 pb"^ of data collected with the 
LHCb detector. The inclusive cross-sections have been measured separately for and 
production as well as the ratios o-w+^^+u/o'w-^fj.-p and {o-w+^^J.+u+crw-^t^-p)/o^z-^^J.^J. 
and the lepton charge asymmetry {aw+^^l+u - aw-^i^-p) /{<^w+^n+u + aw-^ti-p)- The 
results have been compared to five next-to-next-to-leading order QCD predictions 
with different sets for the parton density functions of the proton and to one cal- 
culation at next-to-lcading order. There is general agreement with the predictions, 
though some of the PDF sets overestimate the ratios of the cross-sections. The ratio 
aw+^n+u/cyw-^f^-u =1.27 ± 0.02 ± 0.01 is measured precisely and allows the Standard 
Model prediction to be tested with an accuracy of about 1.7%, comparable to the uncer- 
tainty on the theory prediction. These represent the first measurements of the W and Z 
production cross-sections and ratios in the forward region at the LHC, and will provide 
valuable input to the knowledge of the parton density functions of the proton. The uncer- 
tainty on the cross-section measurements is dominated by systematic uncertainties. Since 
most of these are statistical in nature, the accuracy on the measurement with further data 
is expected to significantly improve. 
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Appendix 

A Tables of results 



Table 2: Correlation coefficients between W and Z in the five bins considered. The 
luminosity uncertainty is not included. 
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Table 3: Differential Z ^ jifi cross-section, daz^nn/dy^ , in bins of boson rapidity. The 
first cross-section uncertainty is statistical, the second systematic, and the third due to 
the uncertainty on the luminosity determination. The correction factor f^sR which is used 
to correct for FSR is listed separately. 
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Table 4: Differential W ^ /lu cross-section, daw-^^iu/v^, in bins of lepton pseudorapidity. 
The first cross-section uncertainty is statistical, the second systematic, and the third due 
to the uncertainty on the luminosity determination. The correction factor which is 
used to correct for FSR is listed separately. 
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Table 5: Lcpton charge asymmetry, A^, in bins of muon pseudorapidity for a pt^ threshold 
at 20, 25 and 30 GeV/c. The first uncertainty is statistical and the second systematic. 
The effect of FSR is at the level of 10~^ and is not listed. 
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Table 6: W cross-section ratio, Rw = (yw^-^ii-^vl <^w-^ii-v-, in bins of muon pseudorapidity. 
The first error is statistical and the second systematic. The effect of FSR is at the level 
of 10~^ and is not hsted. 
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